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Abstract 

A procedure for the construction and the classification of multi- 
lattices in arbitrary dimension is proposed. The algorithm allows to 
determine explicitly the location of the points of a multilattice given its 
space group, and to determine whether two multilattices are arithmeti- 
cally equivalent. The algorithm is based on ideas from integer matrix 
theory, in particular the reduction to the Smith normal form. Among 
the applications of this procedure is a software package that allows the 
classification of complex crystalline structures and the determination 
of their space groups. Also, it can be used to determine the symmetry 
of regular systems of points in high dimension, with applications to 
the study of quasicrystals and sets of points with noncrystallographic 
symmetry in low dimension, such as viral capsid structures. 

1 Introduction 

A monoatomic (A^ + l)-lattice is a regular set of points in M" that is the 
superposition of + 1 identical Bravais lattices, and can be described by a 
reference (or skeletal) Bravais lattice and A^ shift vectors that measure the 
location of A^ points relative to a fixed point of the reference lattice (cf. e.g., 
IMJ). 

The symmetry of a multilattice is determined by those symmetry opera- 
tions of the point group of the skeletal lattice that act as permutations (mod- 
ulo the lattice) on the additional points ([12], [2], [3], [1], [I3], [I5]). More 
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precisely, a symmetry operation of a (A^ + l)-lattice is a triple M = (Mj), 
A = (A^), T = {TD such that, for every a = 1, . . . , iV, 

n N 

Y.M]Pi = Y,PhAi+n, (1) 

0=1 13=1 

with M a unimodular integer matrix which is an element of the lattice group 
of the skeletal lattice, A an integer matrix that describes the permutation 
action of (Mj) on the points of the multilattice, T a matrix of integers repre- 
senting a set of lattice translations, and the components of the a-th shift 
vector. To each triple (M, A, T) a matrix of the form 

can be associated, and it turns out that the set of all triples that satisfy ([1]) 
for a given set of shift vectors is a group under matrix multiplication, the 
symmetry group of the multilattice. 

Denote by r„ the group of all matrices of the form ([2]) for arbitrary M 
unimodular integer, A a linear representation of a permutation on the shifts, 
and T an integer matrix: two (A^ + l)-lattices are arithmetically equivalent 
if their symmetry groups are conjugated in tv- This notion of equiva- 
lence generalizes to multilattices the usual arithmetic classification of simple 
lattices in Bravais types. 

We refer to ([1]) as the master equation of the multilattice. It can be 
used either to compute the shift vectors P, given the symmetry group, or to 
compute the symmetry group given the skeletal lattice and the shift vectors. 

In this work we describe an algorithm to solve the master equation for any 
given symmetry of the skeletal lattice. Precisely, we assign a subgroup of the 
lattice group of the skeletal lattice, construct a permutation representation of 
this group, and solve the corresponding master equation, letting the integer 
translations be taken care of by simple bookkeeping. The procedure is based 
on ideas from integer matrix diagonalization ( [TB] , [S] , [Z] , [ID] , [H] , [I] ) and 
yields automatically a single representative for each arithmetical equivalence 
class of multilattices. 

Hence, this approach provides a basis for an automatic procedure for 
the classification of multilattices with arbitrary number of points, but is also 
yields a simple method to determine regular sets of points in arbitrary dimen- 
sions. This sort of calculations is useful for instance when high-dimensional 
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crystallography is used, via a projection approach, to study quasicrystals or 
sets of points with noncrystallographic symmetry, such as viral capsids that 
display icosahedral symmetry. 

2 Multilattices and the master equation 
2.1 Multilattices 

In this section we introduce the notion of multilattice. We first recall some 
basic definitions about simple lattices, also called Bravais lattices. Unless 
otherwise stated, we use the convention that the sum over repeated index is 
implied. 

Definition 2.1 A Bravais lattice, with basis {ej}j=i...„ (a basis ofW"') and 
origin Qo G M", is the set of points in M" defined by 

C = CiQo,ei) := {Qo + m'Bi G M" / EZfori = l...n}. (3) 

Any other lattice basis {e'-}j=i „ of C is an integer linear combination of the 
vectors ef 

e[ := M^ej M = {Ml) G GL{n, Z) (4) 
with GL{n, Z) the group of integer unimodular matrices. 

Definition 2.2 The point group V of a lattice C is the group of orthogonal 
transformations whose action on a lattice basis corresponds to a change of 
basis of the lattice: 

V={Re 0{n) I 3M = [Mi) G GL(n, Z) / Rta = mIsj} . (5) 

Definition 2.3 The lattice group Q of C is the subgroup of integer unimod- 
ular matrices M that satisfy 

It follows that the matrix representation of V in any basis in M" is conjugated 
to Q in GL{n,R). 

Definition 2.4 Two lattices C and CJ are arithmetically equivalent if the 
associated lattice groups Q and Q' are conjugated in GL{n,Z), i.e., there 
exists H G GL{n, Z) such that: 

g = H-'g'H. (6) 
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Consider now a simple lattice L with basis {ej}j=i^ „ and origin Qo; and 
points Qi, ■ ■ ■ ,Qn not belonging to C and not pairwise equivalent modulo 
L (i.e., such that it does not exist a lattice vector t = m^ei (m* G Z) such 
that = Qp + t ioT a ^ P, and a, /5 = 0, . . . , A^). 

Definition 2.5 A {N + 1) -lattice with basis {ej}j=i^... „ is the union of N + 1 
simple lattices C{Qa, Bi) with origin at the points Qq, . . . , Qat, i-e., the set 

N 

[j{Q^ + m'e,,m' eZ}. (7) 

0=0 

The position of the additional points with respect to the origin of the lattice 
£, called the skeletal lattice, is usually described in terms of the shift vectors 

Pa = Qa~Qo, a = 0,...,N. (8) 

Notice that Pq = 0. 

2.2 The symmetry group of a multilattice 

The symmetric group Sn+i, viewed as the group of permutations on the 
(A^ + 1) points Qo, . . . , Qn, 

a : {Qo, . . . ,Qn) ^ (Qaio), ■ ■ ■ ,Qa{N)), a E Sn+1 (9) 
also acts on the shift vectors in a natural manner. 

Proposition 2.6 Sn+i acts linearly on the 'L-module generated by the shift 
vectors {p^, . . . ,Pn} as follows: 

^ -Pa^ Pa{a)-Pa{0), a = l,...,N, (10) 

for every a G Sn+i- 
Proof. By ([8]) and ([9]), 

Pa = Qa — Qo ^ Qa(a) " Qa(0) = Qa{a) — Qo — {Qa(0) " Qo) 

= Pcr(a) ~ Pcr(0)5 

which proves the thesis. ■ 



4 



We denote by S^^i the group of matrices defined by the hnear action 
(fTUj) of Sn+1 on the shift vectors: 

5;v+i = {{A^J e GL(iV,Z)/3a e ^;v+i : A^^p^ = = - p,(o) } . 

(11) 

The symmetry of a multilattice can be described by the set of triples {R, A^, to), 
with e P, {A^) G iStv+i and ta E C for a = 1, . . . , N , such that the action 
of R on the shift vectors is a permutation of the additional points, modulo 
translations of the lattice: i.e., 

Rp, = Aipf, + t^, a = l,...,N. (12) 

Granted (O, and writing p„ = P^e^ and to, = T^Bi, with e M{nx N, R) 
G (n X N,Z), we may rewrite (fT2l) in the form 

MjP^ = P^A^ + r:, a = l,...,iV, (13) 

i.e., letting M = (Mf ), P = (P^), A = (A^), T = (T^), 

MP = PA + T. (14) 

We refer to (fT3!) or (fT^ as the 'master equation'. We now show that (HM 
defines indeed a group, which completely identifies the symmetry of the mul- 
tilattice. 

Proposition 2.7 Let P = (P^) be given, and let Ti he the subset of Q < 
GL[n, Z) of matrices M = (Mj) such that there exist A = (v4^) G Sn+i and 
T = (T^) E Ai{n X N,Z) that satisfy the master equation [I4\ )- Then 
i) Tt is a group; 

a) the map H S^^i defined by [I4^ is a group morphism that defines a 
permutation representation ofTL on the set {Qo^ . . . ,Qn}, through the rela- 
tion A/^Pp = Q„[a) - Qa{0) ■ 

Proof. Both statements follow from the following argument. If M, H G 
7Y, there exist Am, Ah G ^at+i and Tm,Th integer matrices, such that 

MP = PAm + Tm, hp = PAh + Th. 

^A4{n X N,R) and Ai{n x N,Z) are the vector space and Z-module of 7i x real and 
integer matrices, respectively. 
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Hence, 

{MH)P = M{PAh+Th) = {PAm+Tm)Ah+MTh = P{AMAH) + iTMAH+MTH). 

Further, by multiplying MP = PA^ + Tm to the left by and to the 
right by v4^\ we find 

m-^p = pai}-m-^TmAiI 

and the thesis is proved, since TmAh + MTh and M^^Tm^m matrices 
of integers. ■ 

We denote by r„^Ar the set of matrices in GL{n + A^, Z) defined by 

, , H E 

n,N 



OBJ ^ ' ' V ' ^' ^-j^g^ 

E e M{n X N, Z),B e Sn+i}- 

By the above proposition the symmetry group of a (A^ + l)-lattice can 
then be characterized by the group of matrices JC C r„ jv of the form 

^ = ^'^n,N/Men, MP = PA + tY (16) 

which we shall refer to as the symmetry group of the multilattice. 

Definition 2.8 Two {N + l)-lattices with shift vectors P and P' , and cor- 
responding symmetry groups K, and JC' are arithmetically equivalent if JC and 
IC' are conjugated in the set of matrices Tn,N, i-e-, if there exists a matrix 
Q G Tn.Af such that 

JC' = Q-^JCQ. (17) 

Clearly, Ti is isomorphic to JC but, in general, there are many inequivalent 
multilattices, with different symmetry groups, all isomorphic to the same 
point group H (cf. for instance, [E]). 
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3 The master equation as a system of linear 
equations 

The master equation can be viewed either as a relation that uniquely charac- 
terizes the symmetry group /C of a multilattice, given the shift vectors or 
as an equation in the unknowns that allows to determine all multilattices 
that have a given symmetry group /C. In this paper we take the second point 
of view, and assume that /C, or rather Ti^ is given. Specifically, the problem 
we want to solve is: 

• fix a simple lattice C{Qo, e,) with point group V and lattice group Q; 

• choose the number of additional points {Qi, . . . ,Qn); 

• fix a subgroup 7-^ C ^ of the lattice group of the crystal; 

• choose a permutation representation H Stv+i that associates to each 
M & H a. permutation of the points {Qo, . . . ,Qn)', 

• determine the resulting linear representation 7i Sn+i, in terms of 
the shift vectors; 

• solve the master equation ( fT4l) in the unknowns P for every M G 7Y, 
and corresponding A G ^at+i, and for every possible T. 

• compare the solutions for different choices of T. 
Remarks 

(i) It is enough to look for solutions of (|T3l) such that G [0, 1), i.e., in the 
unit cell of the lattice, defined by 

C(e,) = {p = P'e,/P*G[0,l)} (18) 

for a given lattice basis e^. In fact, given a solution (P*) of (|T3l) . and 5*^ G 
MinxN,Z), then T^^ := T^+Mj5^-5^v4(^ e M{nxN,Z) so that (P^+5^) is 
still a solution of f|T3|) with replaced by T^. The vectors p'^ = p^ + S^^ei are 
translations of the shift vectors p^ by lattice vectors, and thus are equivalent 
modulo the lattice. 
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(ii) Since 7i and /C are finite, they admit a finite set of generators {M'^^\ . . . , M'^^'^) 

and (G(^), . . . , G'(^)), with G^''^ = I ^ ^^^^ j ; Proposition O allows to 

conclude that if the master equation holds for each generator, then it holds 
for all elements of the group /C. Hence ( IT3l) . that should hold for every 
element G of /C, can be replaced by 



k = l K. 



(19) 



3.1 The master equation as a linear system 



Consider first the master equation f[T3|) for a fixed element G 



M T 



A 

/C: it can be rewritten as a conventional system of linear equations. To do 
so, given a G {1, . . . , N} and i G {1, . . . , n}, define 



a = i + (a — l)n. 



(20) 



so that a takes values in {1, . . . , nN}. Conversely, let a = 1, . . . , nN and 
define a and i through the identities 



a 



a — 1 



n 



+ 1, 



a — (a — l)n, 



(21) 



where [■] denotes the integer part of its argument. It is clear that as a 
varies in {1, . . . , nN}, then a and i take values in {1, . . . , N} and {1, . . . , n} 
respectively. Let 

Lt := S^M; - 5]Al (22) 

and 

P'' := P^', := T^, (23) 

where a, i are defined as in (12T|) and, for 6 G {1, ... , nA^} 

^6- r 



/3 



n 



+ 1, 



j = 6-(/3-l)n. 



(24) 



with 5^ and 5j Kronecker deltas. The nA^- dimensional vector [P^) has com- 
ponents that are obtained by ordering of the vectors 

In terms of the vectors P and T and the matrix L, the master equation 
f ll3p takes the simple form 

lapb ^ fa_ ^25) 
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The above assertion follows from a simple argument: let := y|, with 
b = 1, . . . , nN, j = 1, . . . ,n and f3 = 1, . . . ,N consistent with the indexing 
rule above. Then 



nN 



6=1 



_j_ . . . _j_ _j_ yl+(2-l)n _j_ . . . _|_ yn+(2-l)n _|_ 

_^yl+{Af-l)n _^ . . . _^ yn+{Af-l)ri 

Y,' + ■ ■ ■ + Y,^ + Yi + ■ ■ ■ + Y^^ + ■ ■ ■ + + ■ ■ ■ + Y, 

N n 

f3=l j=l 



Hence 



n N N n nN 

E - E ^/^^^ = E E(^i^" + = E ^bP'- 

j=l f3=l f3=l j=l 



b=l 



Consider now the system of master equations (HM for the full set of generators 
of K, i.e., 



k = l,...,K, 



(26) 



with K the number of generators of /C. The associated system of linear 
equations fl2^ is now replaced by a system of the form 



with 

with J = 1, . . . , nNK given by 

J = i + (A; - l)nA^ + (a - l)n, 

with inverse 



(27) 
(28) 

(29) 



k 



a 



J- 1 



+ 1, 



J - (A; - l)niV- 1 



i = J -i^k - l)nN - (a - l)n, 
and the relations between 6, /3 and j are as in ( |2^ . 
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3.2 Tensor form of the master equation 

The relation between the master equation and the matrix L can be rewritten 
in more compact form as follows. For M G GL[n,'L) and A G GL(A^, Z), 
consider the fourth order tensor 

M®/1^, (30) 

with components (M x A^)fp = M'jA^. It is easy to verify that the set of 
tensors of the form (15UI) is a group with product 

{M®A~^)*{H0B~^) = MH(g)A'^B^. (31) 

Further, the indexing rules fl2Ul) and fl^Tl) define a morphism between the 
group of such tensors and the group of invertible nN x nN matrices. 

Proposition 3.1 For i, j = 1, . . . ,n, a, (3 = 1, . . . , N , let 

a = i + (a- l)n, b = j + {/3 - l)n. (32) 
then the association rule 

LI := MiAi (33) 

defines a map M ® AJ ^ L between GL{n, M) GL{N, M), with product *, 
and GL(nA^, M) which is a group morphism. 

Proof. Notice first that if M and A are invertible, then L is invertible, with 
inverse associated to the tensor M^^ ® A^^. Let now := HjB^: then 

nN n N n N 

c=l h=l 7=1 h=l 7=1 

= {MH)i{BAf^ = {MH ® (5A)^);^ = [(M ® A^) * {H ® 5^)]}^, 

(34) 

which proves the assertion. ■ 

Further, the tensors of the form (|30|1 act linearly on the space of real 
matrices Ai{n x N, M) as follows: 

{M ^ A'^) : P ^ MPA, PeM{nxN,R). (35) 

Letting P E M"^ be given by ( !23l) . the above action corresponds to the linear 
action of GL{nN, M) on M"^. In fact, 

nN n N n N 

6=1 j=l 13=1 j=l /3=1 
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4 The solution procedure 



Before discussing (p7|) . we start with a simple remark on diagonal systems. 

Consider first a diagonal system of linear equations with integer coeffi- 
cients 

DX = S, (36) 

with D e M{1 X m,Z) {I > m) and = for J ^ z, X G M*" and S G Z', 
i.e., 

Z)*X* = for 2 < r (no sum on i) 

= for i > r 

with r = rankD. Define 



(37) 



S = {{S\ . . . , S') G Z70 < < for i = 1, . . . , r, = for i = r+1, ...,/} 

(38) 

and 

r 1 2 — 1 1 

0,-,-,...,^^ for. = l,...,r. (39) 



(40) 



Then the solutions X of (l36l) . for 5 G 5, belong to the set 

A" = {(X\...,X™)/X* G forz = l,...,r, 
andX-' G [0, 1) for j = r + 1, . . . , m} . 



Notice that S^^^, . . . , S*' must be zero for the resolubility condition of fl36l) 
to hold, and we have assumed that G [0, 1) for j = r + 1, . . . , m without 
loss of generality. 

The following proposition is a consequence of this simple fact: 

Proposition 4.1 The solutions X of [S^) . with S e 1/ x {0}'"'", belong to 
X mod Z™, I.e., there exists Y e X and K e Z"" such that X = K + Y. 

Proof. Given S' G Z'' x {O}'"*", then for alH = 1, . . . , r there exist Ki e Z 
and Ci E {0, 1, . . . , D'l — 1} such that (no sum on repeated indices) 

S' = D\K, + C\ 

Then L)*X* = D\Ki + (again no sum on repeated indices), and 

X' = K' + Y' with Y' G B\ 
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for i = 1, . . . , r, and the statement is proved. 



Consider now the full system of linear equations fl27j) : instead of solving 
it for a fixed value of the right-hand side ta, and hence of T, we look for 
solutions for every integer vector T. Hence, we rewrite (12 7p in the form 

LiP' = mod Z("^^), (41) 

with J = 1, . . . , nNK and 6 = 1,..., nN and is a matrix with integer 
entries. 

It is a classical result that every such matrix can be reduced to a diagonal 
canonical form, the Smith canonical form , [6j . Precisely, for every matrix 
(lO e M{nNK X nN,Z), there exist matrices [U/) e GL{nNK,Z) and 
{V^'') e GL{nN, Z), such that 

Li = UiDiV,\ (Di) E M{nNK x nN, Z), (42) 

with = for / 7^ a, and D\ divides D\l\ if ^0. The Smith canonical 
form D is unique, whereas the matrices U and V are not. 
Define now: 

V = {P e /P = V~^X with X eX} 
where X is defined as in fl40l) with m = riA^. 



Proposition 4.2 Lei L G A4{nNK xnN,Z) , and D its Smith normal form, 
with r = rank{D): then the solutions of (^I]) belong to V mod Z"^. More 
precisely, there exist exactly oo^"-~^' x [D\D\ . . . D^) solutions modulo 17^^ 
of the system (RIP; and are given by 



P^-V^X'', i^'') - ( J^,J^, ■ ■ ■ ,^,^1, ■ ■ ■ ,tnN-r] , (43) 



with V = V ^ and 

kie{0,l,...,Dl-l}, t,G[0,l). 
Proof. The general procedure to solve fHTj) is as follows: let 

X" = V^P\ 
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so that, since {Uj) G GL{nNK,'L), the system fl41l) can be written in the 
form 

(44) 



DiX" = Q mod Z, 



I.e., 



D\X^ 
DlX^ 








mod Z, 
mod Z, 



(45) 



DrX^ = mod Z, 



where r = rank(Z);f). By Proposition 14. it is sufficient to solve (I45p in the 
set X , i.e. in the set [0, 1)™": we obtain 



— n J- -2_ 

^ ~ '-'5 ni ' ni ' 



X2 = 0, 



^ ^ 



x^'+i = ti e [0,1) 



Di-l 
Dl-l 



(46) 



X"^^ — ^nAT-r G [0, 1), 

with ti real parameters, and this yields ( l43l) . ■ 

Once the X" and the corresponding P'' are computed, the right-hand 
sides of the master equation fl271) may be determined as follows: ffist rewrite 
the diagonalized equation (l44l) in the form 



D^X" 



(47) 



where D^X"^ are known. The solvability condition for (HTll is clearly S""^^ = 
■ ■ ■ = S*"^ = 0, while by construction 5* = fc, for z = 1, . . . , r (cf. Proposition 
14.21) . Once the integer components S"^ are determined, T'^ = UjS^ follows 
immediately. Finally, formulas ([28]) 2,3 yield the solution in terms of the 
and )\ 

The Smith normal form D of a given integral matrix L is unique, but the 
transformation matrices U and V such that L = UDV are not. Since our 
procedure is based on evaluating a fixed solution X of DX = mod Z"^-^^, 
and then reconstructing P as P = V~^X, the question of the independence 
of the solution of U and V arises. We now prove that, for all possible choices 
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of U and V, the solutions belong to the same set, modulo translations of 
the lattice. To see this, recall that, if U* and V* are other transformation 
matrices such that L = U*DV*, then U* = UW and V* = ZV with W and 
Z such that WDZ = D (cf. ^Q\). Let now P* = {V*y^X. Then 

LP* = u*{u*y^L{v*y'^x = u*w-^U'^lv-^z~^x 

= U*W-^DZ-^X = U*DX, 

and since U* is unimodular integer, if DX = mod Z"^^ then U*DX = 
mod JT'^^ . Hence, P* is a solution of the original equation LP* = 
mod Z"^-^.^Then any solution P* = {V*)-^X belongs to V mod Z"^. In- 
fact, since P* is a solution of LP* = mod Z"^^ 



and since U is unimodular integer DZ ^X = mod Z" . Then from 
Proposition 14.11 it follows that Z~^X G X mod Z"^, i.e. there exist G 
Z"^ and X' G A" such that: 

Z-^X = K + X' and P* = V"~^(Z~^X) = 1/"^/^ + 1/-^X', 
so that P* eV mod Z"^. 

5 Arithmetic equivalence 

The following proposition shows that two equivalent multilattices have the 
same Smith normal form, and provides a basis for a criterion to establish 
whether two multilattices are equivalent. 

Consider two equivalent {N + l)-lattices. By definition, their symmetry 
groups )C and /C' are conjugated by some Q G r„ jv of the form 



To simplify things, we assume that the generators (G*-^-*, . . . , G^^^) C fC and 
{G'^^\ . . . ,G''^^'>) C KL' are chosen to be pairwise conjugate, i.e., G'^^'^ = 







LP* = L{V*)-^X = LV^^Z-^X 
UDVV-^Z-^X = UDZ-^X mod Z' 



■nNK 




(48) 
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Q ^G^^^Q for each k = 1, . . . , i^. We write the master equations correspond- 
ing to each multilattice as usual (cf. 



with reduced form 



LP = T, 



DX = 



L'p' = r, 

D'X' = S'. 



Also, for a given square matrix W G M.{nN x nN^ Z), we denote by Wk G 
M. {nNK X nNK, Z) the square matrix of the form 



K 



/ 


w 










w 









\ 









M 









> K times. 



Proposition 5.1 Given two equivalent {N + l)-lattices as above, L and L' 
satisfy the relation ^ ^ 

L' = W^^LW, (49) 



where W G GL{nN, Z) is the integer matrix associated to H ® B ^ in ( [^gp 
through the relation ^33^) . Hence, the matrices L and L' have the same Smith 
normal form: 

D' = D. (50) 

Further, the vectors S, S' G Z"^^ are related through 



S' = U'-^Wr,^US + D'Z, 



(51) 



with U, U' G GL{nNK, Z) are such that L = UDV and L' = U'D'V, and 
Z G Z"^ is a vector of integers. 

Proof. Consider two mutually conjugated generators of K, and /C', 



by hypothesis, 



A^^^ I ' 



A'^'^) I ' 



(52) 
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with Q given by (HHD, so that, in particular, M'^^^ = H'^M^^'^H and = 
B-^A'^^^B. Letting 



then 

In fact, by ([3lD, 

{E'^ ® B~^) * (M(^') ® Jjv - 4 ® (A^'^Y) * {H ® B-^) 

= H~^M^^^H ® B^InB-^ - H-^InH ® B^ {A^^Y 

Assertion fj49l) then follows by letting W be the matrix in GL{nN, Z) asso- 
ciated to W through the rule fl33l) . and using the definition fl28|l of L. 
Further, for each k, ([35]) and ([S2D imply that 

which in turn means that 

f = W^^f + L'J, 

with J e GL{nN, Z) the integer vector associated to the matrix H~^R: J is 
integer by consequence of the fact that H G GL{n, Z) and R e A4(n x N,1,). 
Finally, we obtain flM]) by multiplying the above identity by U'^^. ■ 

When /C has one generator, Wk = in the above proposition, which 
reduces to the statement that for equivalent multilattices L and L' are conju- 
gated. The converse of Proposition 15.11 is not true: as it will be shown later, 
there are multilattices that have the same associated Smith normal form and 
that are not equivalent. 

Proposition 15.11 has a partial converse. 

Proposition 5.2 Given two {N + l)-lattices as above, assume that their 
lattice groups H and Ti' , as well as their images in S^+i, are conjugated, 
i.e., there exists H G GL{n, Z) and B G cSat+i such that, writing 

^ ~ V ac^) j ' ^ ~ ^ J ' 
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for the generators of K, and K,' respectively, then 

^,(fe) ^ H-^M^^^H and A"^^^ = B-^A^^^B, (53) 
for each k = 1, . . . , K . If there exists an integer vector Z G Z"^ such that 

S' = U'-^W^^US + D'Z, (54) 

with the same notations of Proposition \5. the two multilattices are equiv- 
alent. 

Proof. Clearly, fl54|) implies that there exists R E Ai{nx N,1j) such that 

f = W^^f + L'J, 

with J G GL{nN, Z) the integer vector associated to the matrix H^^R which, 
together with (1551) . imphes in turn that 

holds for each k, with Q given by fllS]) . ■ 

6 Applications to the arithmetic classification 
of multilattices 

The above algorithm can help to solve a major problem of the arithmetic 
classification of multilattices, namely to generate all equivalence classes of 
(A^ + l)-lattices. In fact, assume that a subgroup H of the lattice group of a 
Bravais lattice is given. For given A^, there is only a finite number of inequiv- 
alent permutation representations of 7i on + 1 elements. Assume that one 
of these has been chosen, so that a set of generators M'' of Ti, together with 
the corresponding A^''\ are known. This allows to construct the matrix L, 
and Proposition 14.21 yields the shift vectors P and the corresponding trans- 
lations T. To establish which of these solutions are mutually equivalent, we 
can now to apply Proposition [521 with M'(^) = M^'^) and A'^^) = A^''\ In this 
case, the conjugants H and B are just normalizers of (M^^\ . . . , M^^^) and 
of {A^^\ . . . , A^^'^) respectively. To this purpose, it is necessary to compute 
the integer matrices H G GL{n, Z) such that 

H~^M^''^H = M^^\ for all = 1, . . . , iT, 
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and the permutations B G S^^i such that 



for all k 



K. 



We are now in the hypothesis of Proposition 15.21 which yields a necessary 
and sufficient condition for equivalence. 

To illustrate the above procedure, we show that, given a subgroup of 
a lattice group and a permutation representation, our approach yields all 
inequivalent multilattices that have that lattice subgroup and permutation 
representation. Specifically, focusing on 2-lattices, we consider for instance 
the hexagonal point group with trivial permutation representation and show 
how to obtain the inequivalent structures 27 and 28 in [3J . In this case, 
n = ?), N = 1 and K = 2. Consistent with [3], we choose as generators of 
the hexagonal lattice group the matrices 





and A^^^ = A^"^^ = 1. The matrix L corresponding to the master equation is 



with Smith normal form 



D 



I 1 

1 










6 





\ / 
The diagonal system DX 



( -2 
-1 


-2 



V 



-2 
-1 

-2 





3 

-2 
3 





1 
-1 

1 






-2 





1 \ 



-10 

110 
1 






1 / 



V 



1 

3 

-1 




-2 
1 



mod has five nontrivial distinct solutions: 



X, = (0,0,^/6), 
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and the corresponding shift vectors are 



Pi = (2/3, 1/3, 1/2), P2 = (1/3,2/3,0), P3 = (0, 0, 1/2), 
P4 = (2/3,1/3,0), P5 = (1/3,2/3,1/2). 

The translation vectors in the basis that diagonahzes L are 

= (0,0, z, 0,0,0),, ^ = 1,...,5. 

To estabhsh which of these solutions are mutually equivalent, we can now 
apply Proposition [5I2I To this purpose, it is necessary to compute the integer 
matrices H G GL(3, Z) such that 





















1 










1 




\ 





-1 ) 




^ 








and the integers B such that BA^^^B = A^^^ and BA^'^^B = A^'^\ i.e., B = ±1. 
A straightforward calculation shows that H must have the form 

H = ± 

Now, 5*1 and 6*5 satisfy fl5^ . i.e., 

5i = U-^Wj^^US5 + D'Z, 

with B = —1 and H = the identity in M^, which implies that Wk = —Iq 
is the inversion in M^. In fact, for this choice, (IMI) reduces to 





/o\ 




/o\ 




/ 1 


\ 




/ mi \ 

















1 




m2 




1 




5 







6 
















+ 










7714 






















7715 




\^ ) 




\^ ) 






yi 






with irii integers 


, which is satisfied hy Z = 


(0,0,1,0,0,0). 



The same argument shows that 

^2 = U'^W^^USi + D'Z, 
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with the same H, B and Z as before. 

A straightforward check shows that (1541) cannot hold for other choices of 
the normahzers. We conclude that the 2-lattices with shift vectors 

Pi = (2/3, 1/3, 1/2), P4 = (2/3, 1/3,0), 

i.e., the structures 27 and 28 in [3J , are the only inequivalent 2-lattices with 
hexagonal point group and trivial permutation representation. 

In fact, the structure corresponding to the shift P3 is not a 2-lattice: it 
is just the hexagonal Bravais lattice with half vertical lattice parameter. To 
see this, it is enough to notice that, denoting by (ci, 62, 63) the basis of the 
hexagonal Bravais lattice in which the generators of the lattice group 
and M*^^) have the integer representation above, M^^^ and M^"^^ have trivially 
the same representation in the basis (ei, 62, ^63). 
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